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Abstract 

We study the lattice Green's functions of the face-centered cubic lattice 
(fee) in up to six dimensions. We give computer algebra proofs of results 
that were conjectured by Guttmann and Broadhurst for the four- and five- 
dimensional fee lattices. Additionally we derive a differential equation for 
the lattice Green's function of the six-dimensional fee lattice, a result that 
was not believed to be achievable with current computer hardware. 

1 Introduction 

1.1 Bravais Lattices and their Applications 

We consider lattices in M, d that are given as infinite sets of points 

I ^2 nidi : m, . . . , n d G % \ C R d 
^ i=l ' 

for some linearly independent vectors ai, . . . , a<j € R d (throughout this paper, 
vectors are denoted by bold letters). In three dimensions such lattices are called 
Bravais lattices. The simplest instance of such a lattice is obtained by choosing 
a t = Si, the i-th unit vector; the result is the integer lattice 7L d which is also 
called the square lattice (for d = 2), or the cubic lattice (for d = 3), or the 
hypercubic lattice (for d ^ 4). 

The face-centered (hyper-) cubic lattice (fee lattice) is obtained from the 
(hyper-) cubic lattice by adding the center point of each (two-dimensional) face 
to the set of lattice points. In two dimensions this operation is trivial: the 
faces of the square lattice H? are all unit squares with corners (m,n), (m + 
1, n), (m + l,n + 1), (m, n + 1) for integers m, n £ TL. Their center points are 
1? + (i, |) which together with 7L 2 again yields a square lattice, more precisely 
a copy of 1? which is rotated by 45 degrees and shrunk by a factor of The 
situation becomes more interesting in higher dimensions. For example, in three 
dimensions there are 6 faces of the unit cube, and their center points together 
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with all integral translates have to be included. It is not difficult to see that the 
three-dimensional fee lattice consists of four copies of H? , namely 

1? U (1? + (i 1 0)) u {t? + (i o, I)) u {t? + (0, 1 i)) . 

Similarly the fee lattice in four dimensions consists of 7 copies of Z 4 , and in 
general the (i-dimensional fee lattice is composed of 1 + (jf) translated copies of 
TL d . 

The study of Bravais lattices was inspired by crystallography in as much as 
the atomic structure of crystals forms such regular lattices. While the cubic 
lattice is quite rarely found in nature (e.g., in polonium) due to its small atomic 
'packing factor (the proportion of space that is filled when a sphere of maximal 
radius is put on each lattice point, in a way that these spheres do not overlap), 
the fee lattice is more often encountered, for example, in aluminium, copper, sil- 
ver, and gold. The atomic packing factor of the fee lattice is y/2ir/6, the highest 
possible value as was shown by Hales in his famous proof of the Kepler conjec- 
ture [5]. Further applications, in particular in connection with random walks 
(see the next section) appear in physics, chemistry, ecology, economics, and 
computer science when lattice vibration problems (phonons), diffusion models, 
luminescence, Markov processes and other random processes are studied. 

1.2 Random Walks 

For sake of simplicity, the fee lattice as introduced in the previous section, is 
stretched by a factor of 2 in all coordinate directions so that all lattice points 
have integral coordinates. This convention is kept throughout the paper as it 
does not change the relevant quantities that we are interested in (e.g., the return 
probability, see below). 

The aim of this paper is to study random walks on the fee lattice in several 
dimensions. We consider walks that allow only steps to the nearest neighbors of 
a point (with respect to the Euclidean metric). Furthermore it is assumed that 
all steps are taken with the same probability. For example, consider a point 
(k,m,n) in the three-dimensional cubic lattice (2Z) 3 . It is the common corner 
point of 8 cubes. The nearest neighbors in the 3D fee lattice are then the center 
points of some of those faces which have (fc, m, n) as a corner point. Note that 
they all have distance y/2 whereas the other corner and face-center points are 
farther away (their distance is ^ 2) and hence not reachable in a single step. 
Thus the number of possible steps is 8 • 3/2 = 12 (number of adjacent cubes 
times the number of adjacent faces per cube, divided by two since each face 
belongs to two cubes). The same situation is encountered at the center point 
of some face and hence every point in the 3D fee lattice has exactly 12 nearest 
neighbors; this number is called the coordination number of the lattice. 

The above considerations can be generalized to arbitrary dimensions in a 
straight-forward manner; one finds that the set of permitted steps in the d- 
dimcnsional fee lattice is given by 

{( Sl ,..., Sd )e{0,-l,l} d : |«i| + --- + |*«i|=2} (1) 

and thus its coordination number is 4 ( 2 ) . 
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1.3 Lattice Green's Function 



Let p n {x) denote the probability that a random walk which started at the 
origin ends at point x after n steps. Note that in our setting of unrestricted 
walks, c n p n (x) is an integer and gives the total number of walks that end at 
location x after n steps, where c is the coordination number of the lattice. 

In order to achieve information about random walks on the fee lattice, the 
following multivariate generating function is introduced: 

oo 

P{x;z) = Y,Pn{x)z n . (2) 

n=0 

This function is called the lattice Green's function (LGF). By defining the struc- 
ture function A(fe) = A(fci, . . . , kd) of a lattice to be the discrete Fourier trans- 
form 

\{k) = ]T Pi(x)e^ k 

x£TR, d 

of the single-step probability function pi(x), the generating function @ can be 
expressed as the <i-dimcnsional integral 

P(x; z) = — - I ... I — — dfci . . . dk d . 



1 - zA(fc) 



We shall be interested in walks which return to the origin and which we therefore 
call excursions. The LGF for excursions is given by 



„, , , , „ If* f 7 * dfci . . . dkd 

m2) =5> <0)2 -W. -L T^m 



{3! 



In the following, we will only refer to this special instance when talking about 
lattice Green's functions. This function allows one to calculate the return prob- 
ability R, sometimes also referred to as the Polya number, of the lattice. It 
signifies the probability that a random walk that started at the origin will even- 
tually return to the origin. It can be computed via the formula 

Example 1. Consider the square lattice T? which admits the steps (—1,0), 
(1, 0), (0, —1), and (0, 1). Its structure function is 

X(h,k 2 ) = i (e~ lfel +e lkl +e- lk * + e lk *) = 1 (cos k x + cos k 2 ) . 

and therefore its LGF is (see, e.g., [3]) 

1 r f v d/cidfe 2 



P(0, 0; z) = — / / \ ■ * — = -K(^) 

71 JO Jo 1 - | ( cos fc l + cos «2) 1" 

where K(z) is the complete elliptic integral of the first kind. The fact the the 
above integral diverges for z = 1 immediately implies that the return probability 
R = 1; in other words that every random walk will eventually return to the 
origin, a result that was already proven in 1921 by Polya [T5]. □ 
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Example 2. We have already remarked that the two-dimensional fee lattice is 
nothing else but a rotated and stretched version of the square lattice. Neverthe- 
less let's have a look at the LGF when the step set {(—1, —1), (— 1, 1), (1, —1), (1, 1)} 
is taken. Its structure function is 

A(fcl,fc 2 ) = I(V l(fcl+fe2) +e -*( fc l- fc 2) +e i(fcl-fc2) +e *(kl+k 2 )^ 

— — ^ cos(fci + k%) + cos(fci — fe)^ = cos kl COS k2, 

using the well-known angle-sum identity cos(x ± y) = cos & cosy T sin x sin y. 

Although the structure function differs quite a lot from that in Example [T] 
the LGF is the same: 



P(0,0,z) = — / / — = -K{z ) 

K JO J0 — Z COS fel COS K2 7T 



(5) 



as shown in Equation (6) of [4]. Note also that the different distances between 
nearest-neighboring lattice points — 1 in Examplc[T]and s/2 in Example[5] — carry 
no weight since only excursions (and not walks with arbitrary end points) are 
investigated. □ 

It is now an easy exercise to compute the structure function A(fe) for the 
(i-dimensional fee lattice: 

a -l d 
a 



A(fe) = ( ) X] cosfc m cosfc„. 



m— 1 n— m+1 



The LGF is then given as the d-fold integral ([3]) and the return probability can 
be computed by integrating over 1/(1 — A(fc)) 



1.4 The Differential Equation Detour 

The return probability in the fee lattice in three dimensions was first computed 
by Watson [T2] as one of the three integrals which were later named after him, 
and which give the return probabilities in different three-dimensional lattices. 
These probabilities can be expressed in terms of algebraic numbers, 7r, and values 
of the Gamma function at rational numbers. For example, the probability of 
returning to the origin in the 3D fee lattice is given by 

16^4tt 4 

9(r(|)) 6 ' 

For the three-dimensional fee lattice, Joyce [B] expressed the lattice Green's 
function in terms of complete elliptic integrals, too, but the expression is fairly 
complicated and for the higher-dimensionsal fee lattices no such evaluation is 
known at all. Similarly we don't know of any closed-form representation of the 
return probabilities in higher dimensions. 

Instead we will derive differential equations for the corresponding LGFs. 
Although less explicit than the previously mentioned closed-form results, such 
an implicit representation of the LGF provides considerable insight. It allows 
one to efficiently compute the number of excursions for any fixed number of 
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steps, as well as the return probability with very high precision (see Section [3]). 
But also the differential equations themselves reveal very interesting properties 
that are worth investigation. 

To motivate our approach and to illuminate the origin of Equation ((3J) , con- 
sider an arbitrary lattice in TL d with some finite set S C TL d of permitted steps. 
Then clearly the probability function p n (x) satisfies the constant-coefficient re- 
currence 

p„+i(x) = — Vp„(cc-s). (6) 
Let F(y\ z) denote the multivariate generating function 

oo 

Multiplying both sides of (O by y x z n and summing with respect to n and x 
gives 

E E p n +i(x)y x z n = j^E E E^* - 

oo 1 oo 

- z (F(y;z)-l) = ±Y,y SF ^ z ) 



ses 



Thus we obtain 

F(y;z) 



1 ~ M J2 s esy s 

and the lattice Green's function P(0; z) is nothing else but the constant term 
(y°)F(y; z). A differential equation for this expression can be derived from an 
operator 

A(z, D z ) + D yi B 1 + --- + D yd B d (7) 

that annihilates the expression F(y; z)/(yi ■ ■ ■ yd)- Here the symbol D x denotes 
the partial derivative w.r.t. x and the Bj arc differential operators that may 
involve yx, . . . , y^, z as well as D yi , . . . , D Vd , D z . The fact that A may only depend 
on z and D z is crucial and therefore explicitly indicated. In Section [3] we will 
discuss how to find such an operator. From 

(2/1 •••y d )A(g,A) — — + 2^{vi ■■■y d )d V] b =0 

and the fact that the coefficient of y~ 1 in an expression of the form D y E^L-00 a nU r ' 
is always zero, it follows that A((y°)F(y\ z)) = A(P{0;z)) = 0. Also in Sec- 
tion [3] we will see that finding the operator Q is part of computing the d-fold 
integral w.r.t. j/i, . . . , yd- The integral 

dy f f dk 



(yi---yd)(i-^E ses y s ) J J i-zX{k) 

is transformed to © by means of the substitution yj e lkj . 
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2 An Experimental Mathematics Approach 



This section presents some results that were obtained in a non-rigorous way 
using the method of guessing. That is, for finding a linear differential equation 

c m (x)f^(x) + ■■■+ Cl (x)f'(x) + c (x)f(x) = 

for a certain function f(x), one computes the first terms of the Taylor expan- 
sion of fix) and then makes an ansatz with undetermined polynomial coeffi- 
cients Cj[x). If the resulting linear system is ovcrdctermined (i.e., if sufficiently 
many Taylor coefficients were used) but still admits a nontrivial solution, then 
the detected ODE is very likely to be correct. Another strategy to gain confi- 
dence in the result, is to test it with further Taylor coefficients, that were not 
used in the computation. However, this method can never produce a rigorous 
proof of the result and there is always a (very small) probability that the guess 
is indeed wrong. 

2.1 Starting from the Integral Representation 

In this section we briefly recapitulate some previous work done by Broadhurst 
and Guttmann, who used the integral representation ([3]) of the LGF as a starting 
point. 

A differential equation for the LGF of the four-dimensional fee lattice was 
computed by Guttmann [2] (see also TheoremJT]): for this purpose, the four- fold 
integral was rewritten as a double integral whose integrand was expanded as a 
power series. Term-by-term integration yielded a truncated Taylor expansion of 
the LGF which allowed him to apply the method of guessing. 

Very recently, Guttmann [3] reported that David Broadhurst had obtained 
an ODE for the LGF of the five-dimensional fee lattice (see also Theorem [3]), 
a result that required several days of PARI calculations. Broadhurst's strategy 
consisted in expanding the integrand in ([3]) as a geometric scries X)^°=o ^{k) n z n 
and in expanding X(k) n using the multinomial theorem (which gives a (m — I)- 
fold sum if to is the number of summands in A(fc)). The inner terms can now 
be integrated using Wallis' formula 

/W'^=£( 2 ;). 

The structure function of the 5D fee lattice consists of 10 summands. Thus the 
computation of the n-th Taylor coefficient of P(0; z) requires the evaluation of 
a 9- fold sum, or in other words, has complexity 0(n 9 ). 

2.2 Counting the walks 

A different way to crank out as many Taylor coefficients of the LGF as necessary 
is to explicitly count all possible excursions with a certain number of steps. Let 
a n (x) be the number of walks from to a; with n steps and let c denote the 
coordination number of the lattice, then the lattice Green's function P(0; z) = 
Y^ = o a n(0)(z/c) n , as we have already remarked earlier. The values of the 
(d+ I )-dimensional sequence (d again denotes the dimensionality of the lattice) 
can be computed with the recurrence ([6]). To obtain the first n Taylor coefficients 
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hence requires one to fill the {d + l)-dimensional array [a m (x)^ Q<m x x <n 
with values (by symmetry it suffices to consider the first octant only, which again 
by symmetry can be restricted to the wedge x\ xi . . . ^ xa). Still, this has 
complexity 0(n d+1 ). Further optimizations consist in cutting off the regions 
where the sequence can be predicted to be zero (e.g. when Xj > n), and to 
truncate the ccj-coordinates at n/2 (since we are interested in excursions, values 
too far from the origin are not relevant). Although the complexity is better 
than before, computing the full array of values can be quite an effort. Again, in 
the example of the 5D fee lattice, about 115 Taylor coefficients are necessary to 
recover the recurrence for a n (0) (which then gives rise to the differential equation 
of P(0; z)), and hence the full array contains about 115-58 5 /5! rj 6.3- 10 s values! 
Fortunately we can do better. 



2.3 Multi-Step guessing 

How can the recurrence for a n (0) be computed without calculating all the values 
of the multivariate sequence a n (x) in the box [0,n] d+1 (or some slightly opti- 
mized version of it)? In the previous section, we first computed lots of data, then 
threw away most of it, and did a single guessing step. But the guessing can be 
done in several steps which we call multi-step guessing. The method is illustrated 
on the 5D fee example. As before, we start with the recurrence © to crank out 
a moderate number of values for the six-dimensional sequence a n (xi, . . . , xg), 
namely in the box [0, 15] 6 , which takes about 30 seconds only. From this ar- 
ray we pick the values of a n {x\, X2,x 3 , 0, 0) which constitute a four-dimensional 
sequence that we denote with b n (xi, X2, X3). The data is now used to guess 
recurrences for this new function b. One of these recurrences is 

(n + l)b n (xi,x 2 + 3, x 3 + 1) - (n + l)b„(xi,x 2 + l,x 3 + 3)+ 

(n + l)b n (xi + 1, X2, x 3 + 3) - (n + l)b n (xx + T x i + 3, £3)+ 

(n + l)b n (xi + l,x 2 + 3,x 3 + 4) - (n + l)6 n (a;i + l,x 2 + 4, x 3 + 3)— 

(n + l)b n (xi + 3, x 2 ,x 3 + 1) + (n+ l)b n (xi + 3, x 2 + 1,^3)— 

(n + l)6„(xi + 3, x 2 + 1, x 3 + 4) + (n + l)b n (xi + 3, x 2 + 4, a; 3 + 1)+ 

(n + l)b n (xi + 4, x 2 + 1, x 3 + 3) - (11 + l)b n (x 1 + 4, x 2 + 3, x 3 + 1)+ 

(x 2 + 2)b n+ i(xi + l,x 2 + 2,x 3 + 3) - (x 3 + 2)b n+ i(xi + l,x 2 + 3,x 3 + 2)- 

(xi + 2)b n+1 ( Xl + 2, x 2 + 1, x 3 + 3) + (xi + 2)6„+i(xi + 2, x 2 + 3,x 3 + 1)+ 

{x 3 + 2)b n+1 {x 1 +3,.t 2 + l,a;3 + 2) - {x 2 + 2)b n+1 {x 1 +3,x 2 + 2,x 3 + 1) = 

which has the disadvantage that it does not allow us to compute the values 
6„(0, 0, 0) since the leading coefficient vanishes; unfortunately this phenomenon 
occurs frequently in this context. An additional recurrence that does not suffer 
from this handicap is much larger and therefore not reproduced here. Anyway, 
guessing these recurrences can be done in less than a minute. 

Now these recurrences can be used to compute more values for the se- 
quence b n {x\, x 2 , x 3 ) (in 30 seconds one can now go up to n = 30) which in 
turn are used to guess recurrences for b n (xi, X2, 0). These latter recurrences 
allow to compute a„(0) = &„(0,0, 0) for < n ^ 115 in about 2.5 minutes. 
Voila! The whole computation takes less than 5 minutes on a modest laptop. 
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3 A Computer Algebra Approach 



Again, wc want to emphasize that the results presented in the previous section 
are certainly nice, but lack mathematical rigor. To achieve ultimate confidence 
in their correctness we have to apply a different method. One possible such 
method is called creative telescoping, a short introduction of which is given in 
the following section. After that we are able to state our results in the form of 
theorems. 



3.1 Creative Telescoping 

This method has been popularized by Zcilberger in his seminal paper |15j . Since 
then it has been applied to innumerable identities involving hypcrgeomctric 
summations, multisums, integrals of special functions, and various other kinds 
of problems. The basic idea is very simple and we illustrate it on the example of 
a definite integral F{z) = f(x, z) dx. The main step in the algorithm consists 
in finding a partial differential equation for f(x, z) that can be written in the 
form 

(A(z, D z ) + D x B{x, z, D x , D z )) {f(x, z)) = (8) 

where the principal part A € (C(z)(D z ) and the delta part B £ <E(x, z)(D x , D z ) 
are differential operators, with the previously introduced notation of D x being 
the partial derivative w.r.t. x. By <D(z)(D z ) we denote the non-commutative 
Ore algebra that can be viewed as a polynomial ring in the "variable" D z with 
rational function coefficients. The result of the algorithm is a (possibly) in- 
homogeneous linear ODE for the integral F(z) that is obtained by integrating 
Equation ©: 

A(F(z)) + \B(f(x,z))} X ~ b . 

L J x—a 

In applications one frequently encounters the situation that the second part van- 
ishes, yielding a homogeneous ODE. This is because many integrals that occur 
in practice, have natural boundaries. With our study of LGFs, we are in a sim- 
ilar situation: in Section [1.41 it was shown that the principal part automatically 
annihilates the LGF. Nevertheless, in some of the present cases we did do the 
additional (but superfluous) check that the differential equation is indeed ho- 
mogeneous by plugging in the boundaries of the integral, and got confirmation. 

If this method is applied to a one-dimensional integral with hyperexponcntial 
integrand (i.e., its logarithmic derivative is a rational function), then it is called 
the Almkvist-Zeilberger algorithm. Its summation counterpart is the celebrated 
Zcilberger algorithm for hypcrgeomctric summation. For our purposes we have 
to generalize the input class for the integrand to the so-called 9-finite holonomic 
functions: a function f(x\, . . . , Xd) is called 9-finite if for each Xi there exists a 
linear ODE for / with respect to xi. If in addition / is holonomic (the definition 
of this notion is somewhat technical and is omitted here) then the existence of 
creative telescoping operators like ([7]) or ([8]) is guaranteed. 

The first algorithm to compute ([5]) for general 9-finitc functions (our exam- 
ples fall into this class, too) was proposed by Chyzak [T]. It can deal with single 
integrations only and thus has to be applied iteratively for multiple integrals. 
Its main drawback is its complexity that makes it impossible to apply it to the 
problems discussed in this paper. In [9] we have developed a different approach 
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to compute (J5J) which is much better suited for large examples involving multidi- 
mensional integrals. In addition, it can directly deal with multiple integrals by 
computing operators of the form ([7]) , but in the present context it turned out to 
be more efficient to do the integrations step by step. Both algorithms are imple- 
mented in our Mathematica package HolonomicFunctions [lOj . and a detailed 
introduction into the topic is given in [5] . The following example demonstrates 
how this method is applied to the previously discussed two-dimensional lattice. 

Example 3. Looking at the integrand of Equation §5§ one realizes that it is not 
9-finite since no linear ODE with respect to k\ can be found (and analogously 
for fc 2 ). But by means of the simple substitutions cosfci — > x\ and cosfe — » x 2 
we can overcome this trouble: the integral now reads 

P{z) = J_ /' f 1 dx ^ dx2 . (9 ) 

7T 2 7-1 7-1 (1 - ZXlX 2 )\fl - x\^/\ - X\ 

Let f(xi, X2, z) denote the above integrand; it is easily verified that it is a 9-finitc 
function. The three ODEs w.r.t. xi, X2, and z are given by the operators 

Gi = {xix 2 z - l)D z + xix 2 , 

G2 = [x\ — V)(x\X2Z — 1)D X2 + (2x\x\z — x\z — X2), 
Gi = (x\ - l)(xix 2 z — l)D Xl + (2x\x2Z — xi — x 2 z), 

so that Gi(f(xi, X2, z)) = for i = 1,2, 3. In this example it is an easy exercise 
to check that the creative telescoping operator 

, 9 x 9 , 9 X 2 — X^X2 x 2 z — x\z 

z(z 2 - I D + 3z 2 - 1 A + z + D x — ^4 + D X2 - 2 - 10 

X1X2Z — 1 X1X2Z — 1 

annihilates the integrand /. Indeed, it can be written as a linear combination 
g(fl~l) n -w(z 2 + l)-3z 2 + l \^ x 2 , . 

X1X2Z — 1 \x\x2z — \y j \x\x2z — \y 

of the previously computed operators. It follows that the double integral © 
satisfies the ODE 

z{z 2 - l)P"(z) + (3z 2 - l)P'(z) + zP{z) = 

whose solution is the elliptic integral K(z 2 ). 

Alternatively, the two integrations can be performed in two steps (the strat- 
egy that will be applied to the higher-dimensional fee lattices). In the first 
step (integration w.r.t. x\) the following two creative telescoping operators are 
found: 

{x\z 2 - l)D z + x\z + D Xl {x\ - l)x 2 

{x\ - l)(x 2 2 z 2 - 1)D X2 + x 2 (2x 2 2 z 2 - z 2 - 1) + D Xl {x\ - l)(x 2 2 - l)z. 

They certify that the integral J ^ f(x\, X2, z) Ax\ is annihilated by the operators 
(x\z 2 - l)D z + x\z and {x\ - \){x\z 2 - l)D X2 + x 2 {2xlz 2 - z 2 - 1). Next the 
operator 

I 1 i\ n 2 1 10 2 -1 \ 1 \ 1 n X2z(x 2 — 1) 

z[z - l)D z + (3z - 1)D Z + z - D X2 2 2 — — 

(x 2 z ij 

which is a linear combination of the previous ones, again reveals the same ODE 
for the double integral. □ 
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3.2 Results 



Using the above methodology and software, we have computed differential equa- 
tions for the LGFs of the fee lattices in four, five, and six dimensions, and rig- 
orously proved their correctness. Additionally, this allows the computation of 
the return probabilities in the respective lattices up to very high precision. 

Theorem 1. The lattice Green's function of the four- dimensional face- centered 
cubic lattice 



-i P"K /*7T fH fTT 



dki d&2 d&3 d&4 



| ( cos k\ cos &2 + cos k\ cos fc 3 + ■ ■ • + cos k 3 cos k^j 
satisfies the differential equation 
(z - l)(z + 2)(z + 3)0 + 6)(z + 8)(3z + 4) 2 z 3 p( 4 ' (z)+ 

2(3z + 4)(21z 6 + 356z 5 + 2079z 4 + 4920z 3 + 3676z 2 - 2304z - 3456)z 2 p( 3 ) (z)+ 
6(81z 7 + 1286z 6 + 7432z 5 + 19898z 4 + 25286z 3 + 11080z 2 - 5248z - 5376)zP"(z)- 
12(45z 7 + 604z 6 + 2939z 5 + 6734z 4 + 7633z 3 + 3716z 2 + 224z - 384)P'(z)+ 
12(9z 5 + 98z 4 + 382z 3 + 702z 2 + 632z + 256)zP(z) = 0. 

Proof. Here we give only an outline of the proof. The calculations in extenso 
are provided as a Mathcmatica notebook in the electronic supplementary mate- 
rial [TT] (to be downloaded from http://www.risc.jku.at/people/ckoutsch/fcc/). 

The substitutions cos kj — > Xj transform the integrand of the four-fold inte- 
gral to 

f{x!,...,X 4 ,z) = (11) 

(i - § (xix 2 + xi.t 3 h — h 2:33:4)) • rij=i J 1 - 

This expression is 9-finite and thus a Grobner basis of the zero-dimensional 
annihilating left ideal can be computed (annO in the notebook). Next, some op- 
erators A ] (x 2 ,X3,X4,z 1 D X2 , D X3 ,D Xi ,D z ) and Bj ( X\ , ■ ■ ■ , X4, z, D Xl ,...,D Xi ,D z ) 
for 1 < j 4 arc computed, such that Aj + D Xl Bj is an element in the left 
ideal generated by annO. This fact can be easily tested by reducing it with the 
Grobner basis: the remainder being answers the membership question in an 
affirmative way. In the notebook, the A,'s are collected in the variable annO, 
and the P/s in the variable deltal. We conclude that A\, A2, A3, and A4 
generate an annihilating left ideal for the integral J Q f{x\, x%, X3, 2:4, z) dxi. In 
a similar fashion, the integrations with respect to 2:2, 2:3, and 2:4 are performed, 
yielding a single ODE in z that annihilates P(z). □ 

Note that this theorem confirms the result given in [2]. Guttmann also ob- 
served that the differential equation given in Theorem [T] has maximal unipotent 
monodromy (MUM), i.e., its indicial equation is of the form A" and hence has 
only as a root, and additionally satisfies the Calabi-Yau condition. Many 
LGFs of other lattices fall into this class, too, and therefore this fact may not 
seem too surprising. 

In order to receive the return probability in the four-dimensional fee lattice, 
holonomic closure properties are applied to compute a differential equation for 



P(z) = 

n=0 \k=0 
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which in turn gives a recurrence for f(n) = ^2k =0 Pk(0)'- 
(n + 2)(n + 3) 2 (n + 4)(35n 2 + 420n + 1252)/(n)+ 

(n + 3)(n + 4)(595ti 4 + 11375n 3 + 79874?! 2 + 244384n + 276024)/(n + 1) 
3(n + 4)(1015n 5 + 24780n 4 + 240253n 3 + 1156976n 2 + 

2769392?i + 2638272)/(?i + 2) + 
(3325?i 6 + 107100n 5 + 1427695n 4 + 10080600n 3 + 39767416n 2 + 

83134488n + 71984160)/(?i + 3)- 
4(2065n 6 + 62580n 5 + 788848n 4 + 5295615?! 3 

40139838?! + 33590844)/(n + 4)- 
12(735n 6 + 25200n 5 + 359282?! 4 + 2725632n 3 

26259960n + 24690708)/(n + 5)+ 
288(35n 2 + 350?i + 867)(n + 6) 4 /(?! + 6) = 0. 

The initial values 



19973086?! 2 - 
11601091?! 2 - 



/(0) = 1, /(1) = 1, /(2) = §, /(3) = g, /(4) 



1637 
1536' 



/(5) = ^ 
M ' 512 



are easily (of course, not by hand!) computed by counting the number of ex- 
cursions of length up to 5. For the return probability 

fl=l-(lim fin))' 1 

we need to evaluate the limit of the sequence f(n). This can be done very 
accurately when knowing the asymptotics of the sequence. We apply the method 
of Wimp and Zeilberger [14] , which has been implemented in Mathematica by 
Kauers [7], and obtain the following basis of asymptotic solutions: 



si(n 
s 2 {n 
s 3 (n 

s 6 (n 



TV 
1 

ri 
1 
n 
1. 



1 

2 
1 
3 
1 
6 
1 
8 

1 - ■ 



1 



1 - 

7 

9^ 



b 

6ri 
5 

2ri 
45 
14n 
52 
9^ " 



67 



1459 



24n 2 144n 3 



01 i 



2?! + 8n 2 8n 3 Vn 4 //' 



4633 112407 



392?! 2 5488n 3 
812 45820 



O 



27?i 2 



7 



18n 3 



O 



243n 3 
1 



Obviously, the first four solutions do not significantly contribute as they tend 
to very rapidly. Thus by taking into account Ss(n) and se(n) only, and com- 
puting the asymptotic expansion to a higher order (e.g., 30), allows us to obtain 
(at least) 100 correct digits of the limit. 

Corollary 2. The LGF of the four- dimensional fee lattice, evaluated at z 
is 



1 



P(l) w 1.10584379792120476018299547088585107443954623663875285836499, 
and therefore the return probability is 

0.09571315417256289673531676490121018570070881963801735768774. 
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Theorem 3. The lattice Green's function of the five- dimensional face- centered 
cubic lattice 



satisfies the differential equation 

16(z - 5)0 - l)(z + 5) 2 (z + 10)0 + 15 )( 32 + 5)(15678z 6 + 144776z 5 + 
449735z 4 + 933650z 3 - 1053375z 2 + 3465000z - 675000)z 4 p( 6 )(z)+ 
8(z + 5)(3057210z 12 + 97471734Z 11 + 1048560285z 10 + 3939663705z 9 - 
4878146975z 8 - 87265479875z 7 - 304623830625z 6 - 266627903125z 5 + 
254876515625z 4 - 1289447109375z 3 - 503550000000z 2 + 
1774828125000z - 354375000000)z 3 p( 5 > (z)+ 

10(27279720z 13 + 923795772z 12 + 11725276842Z 11 + 68439921540z 10 + 
148313757125z 9 - 382134335775z 8 - 3351125770500z 7 - 7801785421250z 6 - 
3779011321875z 5 - 7716298734375z 4 - 39702348750000z 3 + 
3393646875000z 2 + 23905125000000z - 5568750000000)z 2 pW(z)+ 
5(255864960z 13 + 7892060544z 12 + 92744995638Z 11 + 524857986060z 10 + 
1350059072325z 9 - 465440555100z 8 - 13545524756500z 7 - 26918293320000z 6 - 
3649915059375z 5 - 77498059625000z 4 - 190176960000000z 3 + 
40530375000000z 2 + 45343125000000z - 13162500000000)zP( 3 >(z)+ 
5(496679040z 13 + 13819981248z 12 + 149186684934Z 11 + 810956145330z lo + 
2287368823475z 9 + 1646226060075z 8 - 8282515456375z 7 - 6199228765625z 6 + 
13367806743750z 5 - 110925736437500z 4 - 133825053750000z 3 + 
44457862500000z 2 + 5055750000000z - 3240000000000)P"(z) + 
10(167064768z 12 + 4143853440Z 11 + 40678130502z 10 + 209673119160z 9 + 
607021304825z 8 + 689643286650z 7 - 135661728250z 6 + 3711617481250z 5 + 
2664478321875z 4 - 21210430812500z 3 - 7268326875000z 2 + 
4816462500000z - 189000000000)P'(z)+ 

30(7525440z n + 163913184z 10 + 1443544710z 9 + 6925739310z 8 + 
19123388575z 7 + 21336230625z 6 + 36477006875z 5 + 187923165625z 4 - 
55567000000z 3 - 346865625000z 2 + 84037500000z + 27000000000)P(z) 

Proof. The proof is very analogous to that of Theorem [T] and is given in detail 
in the supplementary material [11| . □ 

Remarkably enough, the indicial equation of the differential equation pre- 
sented in Theorem[3]is A 5 (A — 1) and hence the ODE lacks MUM. For the same 
reason it is not a Calabi-Yau differential equation. 

Corollary 4. The LGF of the five-dimensional fee lattice, evaluated at z = 1 
is 

P(l) w 1.04885235135491485162956376369999275945402550465206640313845, 

and therefore the return probability is 

P« 0.04657695746384802419337442059480329107640239774632112930532. 

Theorem 5. The lattice Green's function of the six- dimensional face-centered 
cubic lattice 



W-i /*.../ 

7T 5 JO JO 



dki d/?2 d&3 d&4 dk§ 



1 — ( cos ki cos ki + cos k\ cos fc 3 + ■ ■ • + cos k^ cos k^) 




1 — ( cos ki cos &2 + cos ki cos k$ + ■ ■ ■ + cos &5 cos ke) 



dfci dfc2 d&3 dfc4 dk§ dk$ 



12 



satisfies a differential equation of order 8 and with polynomials coefficients of 
degree 43. Its leading coefficient is 

z 6 (z - 3)0 - l)(z + 4)(z + 5)(z + 9)(z + 15) 2 (z + 24)(2z + 3)(2z + 15) 
x(4z+ 15)(7z + 60)poly(z) 

where poly(z) stands for a certain irreducible polynomial of degree 25, and its 
indicial equation is A 6 (A — l) 2 . The full equation is too long to be printed here, 
but can be found in JZ2]/. 

Proof. The proof is very analogous to that of Theorem [T] and is given in detail 
in the supplementary material [llj . □ 

As in the five-dimensional fee lattice, the differential equation of Theorem [5] 
lacks MUM and therefore is not Calabi-Yau. 

Corollary 6. The LGF of the six- dimensional fee lattice, evaluated at z = 1 is 
P(l) « 1.02774910062749883985936367927396850209243990900114872425172, 
and therefore the return probability is 
i?« 0.02699987828795612426936417542619638021612262676239501413843. 

4 Outlook and Acknowledgements 

While the calculations for Theorem [T] and Theorem [3] arc done in a few min- 
utes respectively hours, it was a major effort of several days to compute the 
certificates that prove Theorem [5j they are several hundred McgaBytes in size. 
With the methods described in this paper and with the current hardware, it is 
completely out of the question to attack the fee lattice in seven dimensions. An 
interesting question is whether the pattern that showed up in dimensions four 
to six continues. This would suggest a differential equation of order 10 with 
indicial equation A 7 (A — l) 3 . But who knows? 

For the three corollaries we computed the approximations for the return 
probabilities with more than one hundred valid digits. But we have no clue 
what their exact values are. Cyril Banderier evaluated these numbers up to 
several thousand digits, but also he was unable to identify the closed forms. So 
we leave these questions open, as a challenge for future research. 

Anthony Guttmann aroused my interest in computing LGFs during his talk 
at the 65th Seminaire Lotharingien de Combinatoire in Strobl, Austria. I would 
like to thank him for his encouragement and advice. I am grateful to my col- 
league Manuel Kauers for interesting discussions and for pointing me to the 
trick of multi-step guessing. 
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